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We study a bilayer two-dimension-electron-gas (2DEG) adjacent to a type-II superconductor thin 
film with a pinned vortex lattice. We find that with increasing interlayer tunneling, the system of 
half filling presents three phases: gapped phase-I (topological insulator), gap less critical phase-II 
(metal), and gapped phase-III (band insulator). The Hall conductance for phase-I/III is 2/0 e^//i, 
and has non-quantized values in phase-II. The excitation (response to topological defect, a local 
vortex defect) in these three phases shows different behaviors due to the topological property of the 
system, including fractional charge e/2 for each layer in phase-I. While in the case of quarter filling, 
the system undergoes a quantum phase transition from metallic phase to topological insulator phase 
(with excitation of fractional charge e/4). 

PACS numbers: 73.43.-f, 71.10.Pm, 74.25. Qt 



Introduction Topological excitation and fractional 
charge have attracted considerable attention for decades. 
Jackiw and Rebbi first found topological excitation with 
fractionalized charge in one-dimensional (ID) system of 
a Dirac fermion field coupling with a topologically non- 
trivial Bose field Su et al. provided a nice intuitive 
picture for charge e/2 soliton in polyacetylene chain 
Goldstone and Wilczek discussed the possibility of irra- 
tional charge in 1+lD systems Q. The two-dimensional 
(2D) fractional quantum Hall (QH) system with strong 
correlation and time-reversal symmetry (TRS) broken, 
firstly studied by Tsui et al supports elementary ex- 
citations of the many-body ground state with fractional 
charge and fractional statistics Q. Recently, Hou et al. 
have studied a model of graphene-like structures with 
a vortex configuration in the Kekule modulations of the 
hopping amplitudes. They have found a fractional charge 
e/2 bonded to the vortex without breaking TRS. Similar 
phenomenon was found in a model of 2D square lattice 
with the analogical modulations of the hopping ampli- 
tudes 0] ■ These results indicate that the TRS and strong 
correlation are not the necessary conditions for charge 
fractionalization in 2D. 

Weeks et al. fs!] proposed another experimentally ac- 
cessible 2D weakly interacting system of two-dimensional 
electron gas (2DEG) in the integer QH state adjacent to 
a film of type-H superconductor supplying the quantiza- 
tion of flux in units of i$o {^o — h/e). They found the 
excitations with fractional charge and anyonic statistics 
which can be described by a wave function composed by 
a set of filled one-particle states. Moreover, the system 
could be fabricated in the laboratory . 

As is known, bilayer system may show new interesting 
features compared with those of the corresponding mono- 
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layer system, for example, the bilayer QH system versus 
the monolayer QH system. In this paper, we study the 
bilayer 2DEG adjacent to a film of type-H superconduc- 
tor. The feasibility of the similar system with mono- 
layer 2DEG was argued by Seradjeh et al. in detail j^. 
Our studies show that the system undergoes quantum 
phase transitions from gapped phase-I/ topological in- 
sulator (TI) to gapless critical phase-II/metal, then to 
gapped phase-III/band insulator (BI) with increasing in- 
terlayer tunneling at half filling (here, half filling means 
one electron per magnetic cell; note that a magnetic cell 
is twice of a monolayer crystalline cell). These quan- 
tum phase transitions are not related to any symmetry 
breaking. The difference between TI phase and BI phase 
is that the energy spectrum of the edged system in the 
former has gapless topological edge state. Furthermore, 
we find that the three phases have different Hall conduc- 
tances and excitations as response to the flux defect, i.e., 
an extra or missing i<I>o flux in the vortex lattices shows 
different charge density proflles. The excitation of the 
BI phase has charge e/2 for each layer. While at quarter 
filling case, there are one gapless metallic phase and one 
gapped TI phase with excitation of charge e/4 for each 
layer. 

Lattice model and energy spectrum The tight- 
binding Hamiltonian for independent electrons in the 
presence of square vortex lattice is given by 

Hlatt = ~ tije '''"cf^Cja — tiije '''"^"C^Cjo; 

<i,]>a <<i,j»a 

where a{(3)—\,2 refers to different layers, <i,j> 
{<<i, j>>) represents the nearest-neighbor (next- 
nearest- neighbor) sites, tij (tiij) is hopping amplitude 
between nearest-neighbor (next-nearest-neighbor) sites 
in the same layer, and t2 is the interlayer tunneling am- 



plitude. In this work, for simplicity, we assume that tij=t 
and tiij=ti. We have neglected the electron-electron in- 
teraction and also assumed that the spins of all electrons 
are polarized along the field. Cia annihilates an electron 
at site in layer a. In the following we set the lattice 
constant a=l, t=l, t\/t=^i, and ^2/^=72- 

The effect of magnetic field is included through the 
Peierls phase factors 



27r 
$0 



— / A„ • dr. 



_ 



Aq • dr. 



<i,j>a 



«i,3» 



with AQ,=^(0,a;) the vector potential in Landau 
gauge so that each plaquette in layer a is uni- 
formly threaded by a flux ^$o- The magnetic 
unit cell in layer a includes two sites {I, m)a and 
{l+l,m)a denoted by Aa and Ba. In bilayer sys- 
tem, the wave function is of the form of spinor field 

tp{r) = {cBl{ri + l,n),CAlirim),CB2iri + l,n),CA2{rhn), V- 

The bilayer Hamiltonian can be diagonalized in the 
momentum space, with the reduced Brillouin zone 
BZ={k : |fc,|<V2,|fcJ<^}, as Hiau=E).^4f^k^k, 
where Vk=(cBi(k), CAi(k), CB2(k), CA2(k))'^. Further- 
more, after a rotation of the spinor field, </'k='5'~''V'k, 
with 5*= exp(i7r/®(T2:/4) exp{iTrI^az/4:) exp{iTrI^ax/2), 
with (Tx.y,z the Pauli matrices, the Hamiltonian _ffk has 
a simplified form 

iJk = S^HS = / (g) (2 cos kyffx + 2 cos k^cTy 

+ 471 sin fc^ sinfcycr^) -f 72CT2: (g) /. (2) 

We define the operator G=i(jy®I and it is easy to see 
that GH*G=H, G^= —1. As a consequence, for each 
eigenstate ipE with eigenvalue E, there is a corresponding 
eigenstate G^p'^ with eigenvalue —E. The energy bands 
of -ffk are given explicitly by 



£^i(k) =72+e(k),^2(k) =72-e(k), 
EsO^) = -72 + £(k), S4(k) = -72 - e(k). 



(3) 



where £(k)= 2 y cos^ kx + cos^ ky + A-yf sin^ kx sin^ ky. 
The energy bands are symmetric about zero energy. 
The bands are sketched in Fig. 1 for different values 
of interlayer tunneling 72. There are four branches of 
the curves corresponding to Eq. (4). For the half filling 
case with 0<72<47i, the bands i?2(k) and -E3(k) are 
separated with a gap A=87i— 272. The system is in a TI 
phase due to the existence of gapless edge states (with 
topological winding numbers) for the edged sample 
(the detailed calculation is not shown in this paper), or 
nontrivial Chern number as seen later. When 47i<72< 
2V2 [H, the bands E2{k) and Esik) mix each other, 
leading to disappearance of the gap. For 72>2-\/2, the 
gap reappear, showing the behavior of a usual BI. For 
the case of quarter filling, when 72 < \/2—2yi, the bands 
E3 and E4 mix each other leading to gapless metallic 
phase; when 72 > -^2— 271, the system is always in the 
TI phase with a gap A=2(72-(y2-27i)). 
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FIG. 1: (Color online) The energy band structures of the 
lattice Hamiltonian (1) with different 72 : (a) 72 = 0.2, (c) 
72 = 1.0, (e) 72 — 3.0. 71 — 0.1. The corresponding densities 
of the states in arbitrary unit are shown in (b), (d) and (f). 



Quantum phase transitions We look at the quan- 
tum phase transitions with changing of the interlayer tun- 
neling from further other point of view. First, we rewrite 
the system's Hamiltonian with 72 as the control param- 
eter, 



ffk(72)-<' 



-72i/^^) 



(4) 



where H^^=I(g,{2 cos kyax+'2. cos kxay+iji sin kx sin kyaz) 



5/. Note that [H^\h^^]=0. Thus quan- 



and h'^^—Gx 

turn phase transition may appear as a few lowest energy 
levels cross and the properties of ground state change 
dramatically [l^ . One may also find possible quantum 
phase transitions by calculating the ground-state fidelity, 
defined as the overlap between ^'0(72) and \E'o(72+'5) 



f^(72,<5) = 1(^0(72 + <5)|*o(72))|, 



(5) 



where 5 is a small c^uantity, and the many-body ground- 
state wave functions of the system can be constructed 
with the one-particle wave functions through the Slater 
determinant. In the numerical calculation of the fidelity, 
we have used 40000 basis and (5=0.005. The fidelity as 
a function of 72 is shown in Fig. 2. We see that it is 
always constant and equal to unity in the two regions 
of 0<72<47i and 72 > 2^/2 for the half filling case, indi- 
cating that each region corresponds to a specific phase. 
The fact that the fidelity is equal to zero in the region 
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47i<72<2\/2 tells us the system is in the critical phase 
in this region, which is consistent with the absence of 
gap and lack of characteristic length scale in phase-II. 
Similarly, in the case of quarter filling (1^=1/4), there 
are only two phases separated by the critical point 72, c= 
\/2— 271. These results are consistent with what we have 
found based on the analysis of the energy spectrums and 
can be described by the energy-level crossing 12]. 

Topological excitations and factional charges 
We first calculate the Hall conductance based on the 
Kubo formula aH=\C with C—i^Y^ J (Pkz ■ (Vk x 

n 

A„(k)), where the upper limit means that the integration 
is over all occupied states below the Fermi energy Ep. 
Here ^„(k)=i (?/;„(k)| VkV'n(k)) is the Berry phase con- 
nection for nth band with wave function '0„(k). When 
the Fermi level is in the gap, one has the well-known 
(Thouless, Kohmoto, Nightingale, and Nijs) TKNN for- 
mula [l3| C^^^Cn, where C„ — ^ j^^cfikz ■ (Vk x FIG. 2: The fidelity versus inter-layer coupfing 72. 71 = 0.1. 
A„(k)) is the Chern number. The results for the Hall 
conductance is shown in Fig. 3, which reflect the differ- 
ent properties of the three (two) phases for the case of 
half (quarter) filling. It is clearly seen that in the half 
filling case, the Hall conductance is 2/0 (in unit of /h) 
in the phase-I/III, and is not quantized in the gapless 
phase-II; in the quarter filling case, the Hall conductance 
is 1 in the gapped phase-II, and is not quantized in the 
gapless phase-I. 

The topological excitation appears as the response to 
the perturbation added to the system. Here the per- 
turbation is in a form of a flux defect yy'&o {v^^^) to 
the vortex lattice as discussed in Ref. [8]. The vec- 
tor potential in each layer has the following modulation 
5Ka—^^^ (r X z). Correspondingly, the modulation to 
Peierls factors is 

_ 27r 




5Aa ■ dr =r] J dO 

<i,j>a <i,j>a 

S(Pija = / ^-^a ■ dr =7] / dd. (6) 

<<i,j>>a <<i,j>>a 



In a stringy gauge, the Peierls factors can be specified as 
dOij = TT {S(pij = tt) if the string, originating from the 
defect and ending at a boundary, cuts the i-j bond, and 
zero otherwise. These extra Peierls factors have the prop- 
erties f SOij=Tr mod 27r S(pij=TT mod 27r) for closed 
loop containing the defect. They look like a vortex pro- 
file. These topologically nontrivial Peierls factors have 
profound effects on the behavior of the excitation. 

Let us look at the basic properties of the excitation, 
as the response of the system to the flux defect, in the 
three phases for the case of half filling. Since phase-I 
and-III are gapped, thus are incompressible, the charge 
for each layer due to the presence of the flux defect is 
SQ=riaHe/2, as argued in [8|. Combining with above re- 
sults of the conductance, we have SQ=e/2 (for 77=1/2) in 
phase-I and SQ=0 in phase-Ill. To confirm this result and 



FIG. 3: The HaU conductance versus inter-layer coupling 72. 
71 =0.1. 



gain more understanding of the properties of the three 
phases, we perform a numerical calculation based on ex- 
act diagonalizations of lattice Hamiltonian. We have ap- 
plied this method to the systems of sizes up to 48 x 48 at 
different values of interlayer tunneling 72 and obtained 
the charge density profile (per layer) as presented in Fig. 
4. The two layers have the same results since the defects 
in each layer are superposable along the magnetic field 
direction. By integration of the density profile, we get the 
charge for each layer 5Q~e/2 (0) for the gapped phase-I 
(HI), agreeing with the above arguments based on the 
incompressibility and Hall conductance. The situation is 
quite different for the gapless phase II. The charge profile 
is nonlocal and there is no well-defined localized charge 
due to the lack of characteristic length scale in the gap- 
less critical phase-II. This is unlike the gapped phase-I, 
where the presence of localized charge is due to the non- 
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FIG. 4: (Color online) The charge densities of the excitation 
for different 72 : (a) 72 = 0.2, (c) 72 = 1.0, (e) 72 = 3.0, in a 
48 X 48 lattice system with periodic boundary condition, (b), 
(d), (f) are the corresponding planforms. In all cases, we have 
subtracted the background charge density which is uniformly 
equal to 0.5 on each site of the lattice at the half filling case. 
71=0.1. 

trivial Peierls phase profile and characteristic length scale 
determined by the gap. 

Having seen the essential difference between the 
gapped and gapless phases, we now explore more 
difference between the two gapped phases. From the 
different values of Hall conductance, we have already 
seen that the difference between phase-I and phase-HI 
is due to their different topological properties. In fact, 
explicit calculation of edged systems (not shown here) 
shows the existence of gapless edge states in phase-I, 
while they are absent in phase-Ill. More understanding 
can be obtained by studying their low energy behavior. 
In phase-I, the low energy excitation can be described 
by the system Hamiltonian expanded around the points 
K={l,±^),Hji=m{2kya^+2ko,ay+ A'yia,)+ j2cra:®I 
with energy spectrum Ek= ±72±2\/fc2-f 47J ^j^g 
momentum relative to K. In this case, the Hamiltonian 
is of the Dirac type with first-order differential operators. 
When the system is in the phase-HI, the low energy 
excitation can be described by the system Hamiltonian 
expanded around the points K =(0,7r) [or (0,0)], 
Hfi=I(»{2{l-kl/2)a^+2{l-kl/2)cry)+j2(Jx<»I with 
energy spectrum Eg_— ±72±2Vl— fc^. Clearly, unlike 
phase-I, the Hamiltonian for phase-Ill is a Schrodinger 



type of second-order differential operators. Also, since 
the 7i term disappears in i/g, the low energy physics 
in phase HI is thus insensitive to the vortex profile and 
looks more like ordinary BI. These discussions reveal 
the essential difference between TI (phase-I) and BI 
(phase-HI). 

Now we briefly discuss the case with quarter filling. In 
this case, as shown in last section there are two phases: 
gapless critcal phase-I and gapped phase-II. In phase- 
I the Hall conductance is not quantized (as shown in 
Fig. 3) and there is no well-defined localized charge. 
While in phase-II, the Hall conductance is quantized 
at le^/h. The charge for each layer can be obtained 
as i5(5=ry(7//e/2=l/2xl/2e=e/4, which is verified by our 
numerical calculation (not show here). The system at 
quarter filling with fractional charge e/4 is an interesting 
example obviously different from many other situations 
with fractional charge e/2. 

Before ending our discussion, we address briefiy the ac- 
cessability of the parameters. There are three important 
length scales in our bilayer system: magnetic length Z^, 
the penetration length of the type-H superconductor Al, 

and the interlayer distance d. Ib=\ As— , ^^'^ nm, 

SO Ib has several nanometers when the magnetic field of 
~10 T. The pinned Abrikosov lattice constant a has the 
same scale oi Ib, which makes it easier to realize $o/2 
per plaquette than the natural solid with lattice constant 
of ~0.1 nm, requiring the magnetic field in order of 10"* 
T. As discussed in the Ref. [8], the vortices are well 
separated and the defects are well localized from the es- 
timation of Al. The interlayer distance d can be tuned 
in experiment and is related to the key parameter ^2 (or 
72) in our bilayer system. Unlike other natural systems 
such as bilayer graphene with 72<<1 [T^ , 72 can have 
values in a very large regime by tuning the width of the 
wells (which form 2DEG) in the growth direction, as well 
as the width and height of the barrier between the wells. 

Conclusion we have investigated the bilayer 2DEG 
adjacent to type-II superconductor thin film. We find 
that the system undergoes quantum phase transitions 
between several different phases through modulating the 
interlayer tunneling: gapped phase/ topological insulator, 
gapless critical phase/metal, and gapped phase/band in- 
sulator. Depending on different (topological) properties 
of each phase, the system shows different Hall conduc- 
tance, accompanying with appearance/disappearance of 
topological excitation with fractional charge of e/2 or 
e/4. 
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